Aging generates regular motions in weakly chaotic systems 
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Using intermittent maps with infinite invariant measures, we investigate the universality of time- 
averaged observables under aging conditions. According to Aaronson-Darling-Kac theorem, in non- 
aged dynamical systems with infinite invariant measures, the distribution of the normalized time 
averages of integrable functions converge to the Mittag-LefHer distribution. This well known theorem 
holds when the start of observations coincides with the start of the dynamical processes. Introducing 
a concept of the aging limit where the aging time t a and the total measurement time t goes to infinity 
while the aging ratio t a /t is a constant, we obtain a novel distributional limit theorem of time- 
averaged observables integrable with respect to the infinite invariant density. Applying the theorem 
to the Lyapunov exponent in intermittent maps, we find that regular motions and a weakly chaotic 
behavior coexist in the aging limit. This mixed type of dynamics is controlled by the aging ratio 
and hence is very different from the usual scenario of regular and chaotic motions in Hamiltonian 
systems. The probability of finding regular motions in non-aged processes is zero, while in the aging 
regime it is finite and it increases when system ages. 

PACS numbers: 05.45.Ac, 05.40.-a, 02.50.Ey, 02.50. Cw 



I. INTRODUCTION 

Aging is a concept describing slow relaxation phenom- 
ena in spin glasses p], interface fluctuations in liquid- 
crystal turbulence [5] , blinking quantum dots [3J 13] and 
transports in cells [5] . In nature dynamical processes may 
start at time — t a long before the actual measurement of 
the process begins at t = 0. In aging systems, statistical 
quantities measured in the time interval [0, t] are cru- 
cially affected by the aging time t a [5J [7J . For example 
the distribution of time averaged mean square displace- 
ment was considered in [7J using stochastic tools rele- 
vant to diffusion of bio- molecules in the cell. In contrast, 
for stationary and ergodic processes, statistical proper- 
ties of time-averaged observables do not depend on the 
aging time t a if the measurement time is large enough 
(t — > oo). Here, we introduce the aging limit where the 
aging time t a and the total measurement time t goes to 
infinity and the ratio T a = t a /t (the aging ratio) is fixed 
as a constant. We proceed to show universal statistical 
properties of time-averaged observables in a class of dy- 
namical systems, extending infinite ergodic theory [5] to 
the aging regime. 

Exponential separation of nearby trajectories, i.e., 
chaos is a feature found in many dynamical systems. In 
many generic Hamiltonian systems, the phase space is 
either chaotic or regular. However, closer look at dy- 
namics many times reveals a mixed phase space. This 
means that parts of the trajectories are rather regular, for 
example Kolmogorov-Arnold-Moser tori in phase space, 
while others are chaotic [5]. Determining which type of 
motion depends ultimately on the choice of initial con- 
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ditions. Here we investigate a completely new mecha- 
nism for dual structure of dynamics. We show how, for a 
class of dynamical systems possessing an infinite invari- 
ant measure, dynamics is generically split into regular 
and weakly chaotic. This splitting is controlled by the 
age of the process, and hence is completely different from 
usual scenario. 

In weakly chaotic systems the separation of nearby tra- 
jectories is sub-exponential JTU]. In many cases such sys- 
tems have an infinite invariant density |llj . i.e., a non- 
normalizable density (see details below). Since chaos is 
a pre-condition of ordinary statistical mechanics, it is 
not surprising that usual statistical concepts like ergod- 
icity and stationarity break down when treating weakly 
chaotic systems. It is known that some dynamical sys- 
tems with infinite invariant measures show an aging be- 
havior [12] . Here we consider a dynamical system whose 
evolution started at time — t a with t a > 0, while an ob- 
servation of the dynamical process starts at time t = 0. 
Within the observation time (0, t) an observer evaluates 
time averages of observation functions in dynamical sys- 
tems, and we are interested in the ergodic properties of 
these time averages. From the viewpoint of dynamical 
systems, aging means that a density at t = strongly de- 
pends on the aging time t a . If a dynamical system has an 
invariant probability measure, an initial density denoted 
by p{x,t a ) converges to an invariant density as t a — > oo, 
indicating that the aging ratio does not affect statistical 
properties of time-average observables. However, in an 
infinite measure system, p(x,t a ) does not converge to an 
invariant measure as t a — > oo. 

Ergodic theory states that the time averages of inte- 
grable functions converge to constant values (ensemble 
averages) for almost all initial conditions if the dynamical 
system has an absolutely continuous invariant probabil- 
ity measure [13] . If an invariant measure cannot be nor- 
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malized (infinite invariant measure) , there does not exist 
an exponent a (0 < a < 1) such that J2k=o f( x k)/n a 
converges to a non-trivial constant and ±oo) for an 
integrable function f(x) with respect to an invariant mea- 
sure (see details below), where {xk}k=o,... is a trajectory 
[8J. That is, time averages of integrable functions do not 
converge to a constant. However, Aaronson's Darling- 
Kac (ADK) theorem [Ml [15] states that time averages of 
integrable functions converge in distribution: there exists 
a sequence \a n \ cx n a such that 




< x 



da{0d£, as n — > oo, 



(1) 

where d a (x) is the Mittag-Leffler density of order a [5J. 
In other words, the time average depends strongly on an 
initial position but the distribution of the time average 
converges to a universal distribution (the Mittag-Lcfilcr 
distribution) for almost all smooth initial densities. Here 
we investigate ergodic theory under the conditions of ag- 
ing. We show large differences from the ADK theorem in 
the aging regime where both t a and t are large, and pro- 
pose a limit theorem describing distributions of properly 
scaled sums J2T=o f( x k)- 

Recently, it was shown that infinite ergodic the- 
ory plays an important role in elucidating an intrin- 
sic randomness of time-averaged observables in dichoto- 
mous processes modeling blinking quantum dots [16] and 
anomalous diffusion [T7J dBJ ■ Since aging appears in infi- 
nite measure dynamical systems, it is an interesting and 
important problem to clarify whether distributional be- 
haviors of time-averaged observables are affected by the 
aging ratio. Here, we provide evidence that the aging 
limit where t — > oo and t a — » oo but their ratio T a = t a /t 
remaining finite, plays a crucial role in characterizing 
behaviors of time-averaged observables. In particular, 
we show that the distribution of time averages of inte- 
grable functions is determined by the aging ratio T a using 
weakly chaotic systems. 



II. AGING DYNAMICAL SYSTEMS 

We consider maps T : [0, 1] — > [0, 1] which satisfy the 
following conditions for some 71 G (0, 1): (i) the restric- 
tions T : (0,71) -> (0, 1) and T : ( 7l , 1) -> (0, 1) are in- 
creasing, onto and C 2 -extensions to the respective closed 
intervals; (ii) T'(x) > 1 on (0, 71] U [71, 1]; T'(0) = 1; (hi) 
T{x) — x is regularly varying at zero with index 1 + 1/ct, 
T(x) — x ~ a x 1+1 / Q (a > 0). These maps are related to 
number theory [15], intermittency [2"0"H2"2"] . and anoma- 
lous diffusions [23TI26] . One of the best known examples 
is the Pomeau-Manneville map [20l [21] : 



xt+i = T(x t ) — x t + x 



l+l/a 



mod 1. 



(2) 




FIG. 1. Schematic view of an aging process in the map (J2| 
(a = 0.8). A trajectory is given by a solid line while renewals, 
Xt > 71, i.e., a(xt) = 1, are depicted as red crosses on time 
axis, t = is the time at which a measurement starts. 



point at x = and hence trajectory is trapped in its 
vicinity escaping slowly and then is reinjected back. The 
map and its extensions have attracted wide interest since 
it exhibits intermittency, power law distributed waiting 
times, weak chaos to name a few novel effects. According 
to Thaler's estimation [37], an invariant density p(x) is 
given by p(x) ~ h{x)x^ 1 / a for x G (0,1], where h(x) is 
a positive bounded continuous function on [0, 1]. Notice 
that when a < 1 this invariant density cannot be nor- 
malized, due to its divergence close to x = 0. Here p{x) 
is defined up to a multiplicative constant, which we later 
specify. In what follows, we consider an infinite measure 
system (0 < a < 1). In aging systems, the systems start 
at time — t a before the measurement is started at t = 0, 
where t a is the aging time (see Fig. 1). We assume that 
initial points x_ ta are uniformly distributed on [0, 1]. The 
density po(x) at t = 0, when the measurement is started, 
is given by the density of Xt a , i.e., p(x,t a ). 



III. RENEWAL PROCESSES 

Renewal processes are point processes where interevent 
times of points are independent and identically dis- 
tributed (IID) random variables [25J. Therefore, renewal 
processes are characterized by the distribution of in- 
terevent times of renewals. Let us consider the following 
observation function 



<j{x) 



(x < 71) 

1 (s>7i), 



(3) 



In what follows, we use the map ([2| for numerical simula- 
tions. This famous map has a marginally unstable fixed 



where 71 < 1 attains T(7i) = 1. We call the state a = 
the laminar phase while a = 1 the chaotic phase. In 
the chaotic phase, trajectories Xt show usual chaotic be- 
havior because of the condition (ii). Moreover, the jump 
transformation T"( x ' on x € [71, 1], which is a transfor- 
mation restricted to [71 , 1] , has an absolutely continuous 
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invariant probability measure whereas the original map 
T does not, where n(x) = minjfc > 1 : T k (x) e [71, 1]} 
[2"T] . With the aid of these chaotic properties in the jump 
transformation, trajectories a t — o(xt) can be regarded 
as a renewal process because the interevent times between 
the events a(xt) — 1 are considered to be IID random 
variables. The probability density function (PDF) of the 
interevent times (or residence times of laminar phase) is 
given by [33] 



V>(r) 



tA c 



T — > OO, 



(4) 



where the exponent of the PDF is controlled by the non- 
linearity of the map in vicinity of the indifferent fixed 
point through the maps parameter a and A is a constant, 
which depends on the map. 

Let N t be the number of renewals in the time interval 
(0,i) for a process which started on t a , and with IID in- 
terevent times according to Eq. Q . In non-aged renewal 
theory t a = 0, the distribution of the number of jumps 
N t /t a obeys the Mittag-Leffler distribution of order a 
(< 1) 29J. Here, we review a derivation of the distribu- 
tion and extend it to the aging regime [30] . Let S n be 
the sum of the interevent times (S n — t± + . . . + t„), then 
we have the following relation 



Pr(N t < n) = Pr(S n > t). 



(5) 



We notice that the interevent time PDF Q belongs to 
the domain of attraction of stable laws [29] . In what 
follows, we use the notation Pr ta= o(-) for non-aged pro- 
cesses, and Pr(-;t a ) for aged processes. By the general- 
ized central limit theorem and setting n = t a x, we have 

Pr (N t /t a <x)= Pr (S„/n 1/a > x' 1/a ) (6) 
t a =o t a =o 



-1/0 



l a (y)dy as t 00, (7) 



where l a (x) is the one-sided stable density with index a, 
which depend on A and its Laplace transform is given by 



/ l a (x)e~ sx dx = exp(-r(l - a)(As) a ). 
Jo 



(8) 



In aging renewal processes, the PDF ipo(T;t a ) of the 
forward recurrence time, that is the time between the 
start of an observation t = and the first renewal (see 
Fig. 1), is different from Q. According to [501 GH], the 
double Laplace transform of i^o{t; t a ), 



i>o(s,s a ) 
is given by 



00 /-OO 



^o(r,i a )e-*— - T 'drdt , (9) 



1po(s, So) 



Furthermore, by Dynkin's limit theorem 
PDF (t a — > 00) reads 



^o(T;£a) 



sin(7ra) 



T a {t a + T)' 



(10) 

the limit 
(11) 



The probability of Nt — is given by ^o(t; t a )dr, 
while, for Nt > 1, the probability Pr(N t < n;t a ) is rep- 
resented by the convolution of i/jo( t ', ta) an d Pr ta= o(N T < 
n - 1) 

Pr(N t /t a <x;t a ) = l-m a (T a ) 
+ [ Pr (N t _ T < xt a - l)Mr;t a )dr, (12) 

where m a (T a ) = ^^{r^tajdr, which is represented by 
the incomplete beta function, 



,„ . sin(7ra) „ 
m a (T a ) = — - — -B 



T 



1 — a, a 



(13) 



The 1 — m a (T a ) describes trajectories with no renewal 
events in (0,t), or in the context of dynamics of maps 
particles which did not escape the vicinity of an unstable 
fixed point in the observation interval. We note that T a is 
the aging ratio, defined in introduction. In the same way 
as a calculation of the probability of N t /t a in non-aged 
renewal processes, we have 

Pv(N t /t a <x;t a ) 

-^l-m a (T a )+ / ip {T-,t a ) I ^ l a (y)dydr, (14) 

for t — y 00 . As a result, the PDF p a (£;T a ) of f = N t /t a 
in the aging limit, t a /t — ¥ T a (t,t a —> 00), is written as 



p a (Z;T a )=6{Z)[l-m a (T a ) 



sin(7ra) 



^0^1+ 1/« 



r(i + 2/) 



This result is consistent with [7] . We note that the distri- 
bution depends strongly on T a even when the total mea- 
surement time goes to infinity. For T a —¥ (non-aging 
limit), p a (£,',T a ) converges to the Mittag-Leffler density 
as expected. 

One can show that the mean (Nt) [7J is 



(N t ) = W+J^B 



1 



1, a 



(16) 



T(a) \T a + l' 

It will soon be usefull to define a normalized variable 



X = N t /(N t ). Using Eq. £5} and C A = t a /(N t ), we 
have 



P a (x;T a ) = -^-p a (x/C A ;T a ). 



(17) 



We note that PDF P a (x', T a ) has the advantage of being 
A independent (PDF p a (x',T a ) depends on A). 
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IV. RESULTS 



A. Distributional limit theorem in the aging limit 



The distributional limit theorem ( 14 ) in aging renewal 



processes implies that the time average of a(xk) con- 
verges in distribution: 



Pr ( t\ E ^ < x ) = P W* Q < x) 

\ k=0 / 



(18) 



Pa (Z;T a )dt. (19) 



Because o~(x) is an integrable function with respect to 
an invariant measure, this distributional limit theorem 
is a generalization of ADK theorem ([!]). By the Hopf's 
ergodic theorem [33], the ratio of the sums of arbitrary 
integrable observation functions f{x) and a(x) converges 
to a constant for almost all initial points: 



SLo/frfc) Jo /jg 
Efe=o a ( x k) Jo adp 



as n 



(20) 



where p is an invariant measure. Therefore, we have the 
following theorem. In the aging limit t a /t — » T a as t a and 
t — > oo, for all integrable functions f{x) with respect to an 
invariant measure \x, the time average of f(x) converges 
in distribution: 



Pr 



{-^ E f( x >) < X j ~* [ P ^ 



(21) 



where Cf — f Q fdp/p([ji, 1]). We note that the distri- 
bution depends on the aging ratio T a and a while we do 
not represent an initial density explicitly in the left-hand 
side in ([21 1. 



B. From Dynkin's limit theorem to evolution of 
density 

Here, we give an explicit representation of an initial 
density in aging processes. In aging renewal processes, 
the probability that there is no renewal until time t 
is given by J. ipo(r;t a )dT. Corresponding probability 
in the map pi) is the probability that trajectories do 
not escape from the interval [0,71), which is given by 
J 7t p(x,t a )dx, where T t ('y t ) = 1 for j t < 1- Using a 

continuous approximation, T'(^ t ) ~ T ^!y > t Z^J t 1 +1 ' > ; near 
x = and T(j t ) = 7t-i, we have 



7t_i - 7t 
It - 7t+i 



- 1 



1 



li 



l/a 



(22) 



It follows that 74 ~ a a t a (the rigorous proof is given in 
[27]). Therefore, we have the following relation: 



p(x,t a )dx ~ C' 4> {T;t a )dT, 



(23) 



where C is a constant independent of time. Differenti- 



ating both sides of (23 1 with respect to t and using (11 1, 
we have 



p(a a /t a ,t a )- 



C 



, sin(7ra) t c 



(24) 



t a+1 " 7T t a {t a +t)' 

As the result, we obtain an initial density in the aging 
process: 



p(x,t a ) ~ C- 



1 



(25) 



1 + t a x 1 / a /a' 

The constant C is the normalization constant which de- 
pends on t a , for t a 3> 1 

fl 1 \ 
C = I / dx 



l + taxVa/c 



sin(7ra) / t a 
ira \ a 



Surprisingly, evolutions of the density are in very good 
agreement with the above estimation even in small num- 
bers of iterations and on whole space [0, 1] (see Fig. 
2). The density cannot converge to an invariant den- 
sity (equilibrium density) and will converge to the delta 
function S(x) as t a goes to infinity. This is a direct evi- 
dence of aging in dynamical systems. We note that the 
scaled density converges to an infinite invariant density 
p(x) [33: 

lim tl- a p(x, t a ) = p(x) ~ sm ( TO ) Q -a x -l/o , x ^ Q) 

(27) 

Moreover, by the change of variable, x — a a y/t™, the 
scaled density, q(y) = a a p(a a y/t";t a )/t", gives a uni- 
versal master curve: 



q(y) = 



sin(7ra) 



1 



1 



y 



l/c 



(28) 



This result is consistent with a rigorous result in general 
intermittent maps by Thaler 35J. Figure 3 shows a con- 
vergence of the scaled density to the master curve. We 
note that the master curve is universal in the sense that 
it does not depends on details of the map except for near 
the fixed point x — 0. 



C. Dynamical instability 

To investigate an effect of the aging on the dynamical 
instability, we consider the Lyapunov exponent. In gen- 
eral, weakly chaotic systems with infinite invariant mea- 
sures have zero Lyapunov exponent [1311 > 36, 37 1. How- 
ever, these dynamical instabilities are known as a subex- 
ponential instability quantified by the generalized Lya- 
punov exponent A a |11[ 137] , which is defined as the av- 
erage of the normalized Lyapunov exponent, A Q = (A Q ), 
where (•) is an average with respect to an initial density 
and 



1 



t-i 

Xa = A 111 To E ln \ T '( X k)l 
t— foo t * — • 

k=0 



(t« = 0). (29) 
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FIG. 2. Initial densities for different aging times t a = 10, 10 
and 10 J (a = 0.8). Symbols with lines are the results of nu- 
merical simulations. Blue dashed curves are the theoretical 
curves without fitting parameters. For all times t a , the den- 
sities are in a nice agreement with the theory. 




y 



FIG. 3. Scaled density (a = 0.6). Dashed line is the master 
curve Eq. ( |28[ ) for a = 0.6. The scaled density approach to 
the master curve as t a — > oo. 



In non-aged systems, A a does not depend on an initial 
density with the aid of ADK theorem. To investigate an 
effect of aging, we consider the generalized Lyapunov ex- 
ponent in the aging limit t a /t -> T a , A a (T a ) = (X a (T a )}, 
where 

1 t_1 

X a (T a )= lim (t a /t^T a ). (30) 



fc=0 



The aging generalized Lyapunov exponent is represented 



A a (T a ) = — / / g(x)p(x,t' + t a )dxdt', 
t Jo Jo 



(31) 



where g{x) — ln|T"(a;)|. Using the limit density (27l, we 
obtain 

A*(T a ) = ( —^- Q a - J g(x)p(x)dx. (32) 



for t ^> 1 and t a 1 with T a fixed. Here p(x) — 
limt a _).oo ta~ a p(x, 0) defines a unique infinite invariant 
density. According to ADK theorem, the non-aging gen- 
eralized Lyapunov exponent (t a = 0) is given by |37j 



i r 1 

A a {T a = 0) = - / 
" Jo 



g{x)p{x)dx. 



In the aging limit, this relation is generalized as 



A Q (T a ) = 



(1 + T a )«-T« 



a 



g(x)p(x)dx. 



(33) 



(34) 



This aging effect on the generalized Lyapunov exponent 
has been confirmed numerically (Fig. [4]). The result 
means that the dynamical instability becomes weak as 
the system ages. 

Using the distributional limit theorem (21), we ob- 



tain the PDF of the normalized Lyapunov exponent 
X = K(T a )/A a (T a ), given by PDF g^, P a ( X ;T a ). As 
shown in Fig. 5, the PDF does not depend on the total 
measurement time t if the aging ratio T a is fixed, and the 
strength of the delta peak at x — is increased as the ag- 
ing ratio T a is made larger. This delta peak corresponds 
to trajectories that do not escape from [0, 71] until time 
t, which weakens the dynamical instability. These trajec- 
tories are regular rather than chaotic (not even weakly 
chaotic). In fact, these trajectories can be treated by 
a continuous approximation. For x t = 0, intermittent 
maps we consider can be written as an ordinary differen- 
tial equation [3TJ |53] : 



dx t 
~dt 



1+1/ a 1 r*j n\ 

a x t [x t = 0). 



The solution is given by 

/ -1 —1 l/a ,\- 

x t = xo{l - a x ' a t) 



(35) 



(36) 



for xt <C 1 or t <C ctx /clq. The difference of nearby 
trajectories Xt and x' t such that xq and x — xq + 5x 
(Sx <C xq) is described as 

\x t - x' t \ = 5x(l + xl /a (l + oT^aot). (37) 

In the aging limit, these regular motions appear even 
when time t goes to infinity. Therefore, in the aging 
dynamical systems, regular motions (RMs) and weak 
chaotic motions (WCMs) intrinsically coexist in the ag- 
ing limit. This implies that when t a » t a large fraction 
of particles are located close to the unstable fixed point 
and hence their motion is regular. 
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FIG. 4. Generalized Lyapunov exponent as a function of the 
aging ratio T a (a — 0.8, 0.5, and 0.2). Different symbols are 
the results of numerical simulations for different a, where the 
total measurement time t is fixed as 10 6 . Curves are theoret- 
ical ones without fitting parameter (A a (T a = 0) is obtained 
numerically). Aging strongly affects the dynamical instability 
when the aging time t a is larger than the measurement time 
t. 



V. DISCUSSION 

We have shown the distribution of time averages of 
integrable functions in the aging limit, which is a gener- 
alization of Aaronson's Darling-Kac theorem. Although 
we use intermittent maps, this generalization will be valid 
for all weakly chaotic maps with infinite invariant mea- 
sures (more precisely, a conservative, ergodic, measure 
preserving transformation) because aging in dynamical 
systems implies that the density does not converges to 
an equilibrium density even when time goes to infinity 
(non-equilibrium non-stationary density). In fact, we 
have confirmed the distributional aging limit theorem in 
the Boole transformation [3]. 



VI. CONCLUSION 

The aging ratio T a plays an important role in charac- 
terizing aging systems. In the aging limit, the distribu- 
tion of time-averaged observables converges to a univer- 
sal distribution, which is determined by the aging ratio 
T a and the exponent a characterizing infinite invariant 
measures in dynamical systems. The mathematical ba- 
sis of this universal distribution is both the generalized 
central limit theorem and Dynkin's theorem for the for- 
ward recurrence time. We have also shown how to use 
the infinite invariant density to calculate statistical aver- 
ages like the measure of separation A Q within the aging 



regime. From Eq. (34) we see an aged dependent pre- 
factor (1 + T a ) a - (Ta) a which multiplies the non-aged 
average. Similar averages hold for other time-averaged 
observables, which are integrable with respect to the in- 
finite invariant density. Thus the infinite invariant den- 
10\ 

T=\ (^10" 

r=i </=io 5 ) 

T=0.l (*=10 6 ) 

r=o.i (/=io 5 ) 



10" 



10" 1 



1(T 



RM 



WCM 



10- 



10"' 




FIG. 5. Probability density function P x (x) °f the normal- 
ized Lyapunov exponent (T„ = 0.1 and 1, t = 10 5 and 
10 6 , and a = 0.5). Symbols with lines are the results 
of numerical simulations. Solid curves are the theoretical 
ones calculated by a numerical integration of the convolu- 
tion in ( |15[ l using the stable density with index 1/2, i.e., 
'1/2 (aO = cexp(— c 2 / (2x))/\/2nx 3 , where c is a scaling pa- 
rameter depending on T a . Probability of regular motions is 
increased when the system ages. Almost all trajectories in the 
left-hand side of the red line exhibit regular motions (RMs) 
while trajectories in the right-hand side are weakly chaotic 
motions (WCMs), where the value of the red line (x = 0.05) 
represents the bin size of the PDF with histogram (the bin 
is needed to graphically represent a delta function). Theory 
is in good agreement with the numerical results including the 
delta distribution. 



sity plays an important role for determination of ergodic 
properties of aging processes. 

We have found that the dynamical instability is clearly 
divided into two different instabilities, i.e., regular mo- 
tions and weakly chaotic motions, in the aging limit. Co- 
existence of regular and chaotic motions is reminiscent of 
generic Hamiltonian systems. However, the meaning of 
the coexistence is completely different. In aging dynami- 
cal systems, the probability of finding regular motions is 
increased according to the aging ratio T a , whereas reg- 
ular and chaotic phase spaces do not depend on T a in 
generic Hamiltonian systems. 
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